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Dimensional crossover of ordering in ferromagnetic films with both periodic and free boundary con- 
ditions is studied for the exactly solvable uniaxial model of classical D-component spin vectors in the 
limit D —> oo. Analytical and numerical solution of the exact equations describing this model shows 
that for lattice dimensionalities d > 4, finite-size corrections to the bulk values of Tc are character- 
ized by the mean-field exponents and anisotropy-dependent amplitudes. For d < 4, the mean-field 
behavior is only realized in the region n^N ^ 1, where Kc is the dimensionlesss inverse transverse 
(with respect to the easy ax;is) bulk correlation length at Tc and A'' is the number of layers in the film. 
In the region HcN <^ 1 and the dimensionality range 3 < d < 4, finite-size corrections are described 
by the universality class of the isotropic D = oo model. For d < 3, magnetic ordering vanishes in 
the isotropic limit, Kc 0, since the film behaves as an object of dimensionality d' = d — 1 < 2 
and long-wavelength fiuctuations destroy the order. Here the suppression of Tc of a film can be 
substantial, depending on the competition between the weakening anisotropy and the increasing film 
thickness. For thick films, Tc becomes small only for very small anisotropy. 

PACS number(s): 64.60. Cn, 75.10.Hk, 75.30.Pd, 75.70.Ak 



I. INTRODUCTION 

Ferromagnetic ordering in the film geometry shows a 
dimensional crossover on the film thickness, which is mea- 
sured by the number of layers N for the lattice spacing 
ao set to unity, and on the closeness to the Curie point 
Tc- For a spatial dimensionality d, thick films, ^ 1, 
show a d-dimensional behavior far enough from Tc. On 
the other hand, in the close vicinity of Tc the film be- 
haves as an object characterized by the reduced dimen- 
sionality d' = d — 1. The latter regime is realized for 

^ N, where is the dimensionless bulk correlation 
length. An important implication of this dimensional 
crossover is that isotropic Hamiltonians cannot describe 
ordering of films in three dimensions. Indeed, since films 
behave two dimensionally at the would be critical point, 
the long-wavelength fluctuations (spin waves, Goldstone 
modes) preclude ordering. Due to the exponential in- 
crease of the correlation length with lowering tempera- 
ture, two-dimensional magnets are extremely sensitive to 
the uniaxial anisotropy which stabilizes order. It is intu- 
itively clear that the stabilizing effect of the anisotropy 
in thick films should be much stronger than that in usual 
two-dimensional systems. 

Finite-size effects and crossovers between different uni- 
versality classes in thin magnetic fikaa_lia,"\£e been ob- 
served in a number of experiments.EJjj'LTulj'Ll The sup- 
pression of Tc and other effects in films of arbitrary 
thickness A^ are frequently addressed with the spatially- 
inhomogeneous mean-field approach of Ref. |^, which, 
naturally, misses the influence of Goldstone modes on 
ordering mentioned above. The general large- A^ asymp- 



totic form of the Tc-shift in thick films is given by 

[Tc{(x^)-Tc{N)]/Tc{cx^)^A/N\ (1.1) 

as was initially established by the high-temperature series 
expansions (HTSE) for the Ising modelBcl For the expo- 
nent A the finite-size scaling theoryll3lld yields A = l/i^b, 
where Vf, is critical index for the bulk correlation length. 
The coefficient A in Eq. (1.1) should depend on the 
anisotropy. For spatial dimensions d < 3, it should 
strongly increase with approaching the isotropic limit, 
so th at fo r any thickness N the right-hand side (rhs) of 
Eq. (1.1) eventually approaches the value of 1. In this 
case of a subs tant ial suppression of Tc, the functional 
form of Eq. (LT) is no longer valid. For very small 



uniaxial anisotropy, the film orders at such low temper- 
atures that all spins along the direction perpendicular 
to the surface are strongly correlated and can be consid- 
ered as a single composite spin. Thus the film is mapped 
on d'-dimensional monolayer with the exchange interac- 
tion NJ. For the model with a uniaxially anisotropic ex- 
change, the anisotropy of the exchange interaction 77' <C 1 
is replaced by drj' / {d — 1), since the composite spins ac- 
quire an "internal" anisotropy. In this case one has 

Tc[d, J, 77', N] = Tc[d - 1, NJ, drj'/id - 1), 1]. (1.2) 

For the model with a single-site anisotropy, the latter, 
relative to the exchange, is not renormalized. It seems 
that no one of conventional theories can reproduce the 
crossover between the two regimes described by Eqs. 
( pni ) and (|r2|). 

The effects of would be Goldstone modes on ordering in 
weakly-anisotropic low-dimensional magnetic structures 



1 



can be conveniently studied within the exactly solvable 
anisotropic spherical model (ASM) which is the limit 
Z? — !■ oo of the classical D-component spin- vector model. 
The Hamiltonian of the latter has the form 



D N 

a=2 / 



(1.3) 



where rtii is the normalized Z?-component vector, |mi| — 
1, and r/ = 1 — 77' < 1 is the dimensionless anisotropy 
factor. In the isotropic case, 77 = 1, the partition func- 
tion of this model coincides in the limit D—r^ 00 (Rcf. 
Ill) with that of the spherical model (SM).y There are, 
however, a number of essential differences between both 
models. In particular, the spherical model yields unphys- 
ical negative coefficient A in Eq. (1.1) for the film with 



free bouncLary conditions (fbc) because of the global spin 
constraint Jlil Improved, versions of the SM without the 
global spin constrainlll3't£l yield positive A, as it should. 
However, within the SM the problem can be only con- 
sidered in four dimensions, since the SM cannot be done 
anisotropic and in three dimensions Tc is zero for all finite 
values of N. 

By contrast, the ASM-describes ordering of bulk sam- 
ples in low dimensions J13 since anisotropy gives rise to 
the gap for spin fluctuations. In Ref. |l^ it was shown 
that even in the isotropic case there are different, longi- 
tudinal and transverse correlation functions (CF) in the 
ASM below Tc, in contrast to the inherently single CF 
in the SM. The ASM correlation functions have a kind 
of spin-wave form, which is, however, valid in the whole 
range of temperatures. This makes the ASM a convenient 
tool for the investigation of classical spin systems at el- 
evated temperatures at the level beyond the mean-field 
approximation (MFA). Whereas the critical coupling of 
fluctuations dies out in the limit D —i- 00 and the criti- 
cal behavior is simplified, the not less important qualita- 
tive effects of the would be Goldstone modes on ordering 
are accounted for in the "pure" form. The most recent 
example of the application of the ASM to the spatially- 
homogeneous systems is Ref. ^ where the spin-liquid 
state of the classical antifcrromagnet on the kagomc lat- 
tice is considered. 

In the inhomogeneous case, the closed lattice-based 
system of equations describing the ASM was obtained 
in Ref. This system of equations has been analyti- 
cally solved in the domain-wall geometry for the biaxial 
generalization of the Hamiltonian (1.3). This model fea- 
tures a phase transition of an interface, which can be 
described analytically beyond the MFA and observed in 
experimented A more complicated situation is realized 
in the semi-infinite geometry, where analytical solutions 
can be only obtained in particular and limiting cases and 
application of numerical methods is necessary. In Ref. 
the correlation functions of the semi-infinite ferromagnet 
have been computed at T >Tc in the case of the so-called 
ordinary phase transition (the transition at the surface, 
that is characterized by its own set of critical indices, 
is driven by the phase transition in the bulk). For the 



film geometry, the problem becomes even more difficult 
because of the additional effect — suppression of Tc in 
comparison to its bulk value. In my preceding Letter 
on ferromagnetic films, Ref. only some analytical re- 
sults for ^ 1 have been obtained, but no numerical 
calculations have been performed. 

The aim of this paper is a complete solution for the 
correlation functions and Curie temperatures of ferro- 
magnetic films within the ASM in the whole range of pa- 
rameters with the use of the numerical algorithm of Ref. 
p^ . In contrast to Ref. ^ calculations will be performed 
for the model with a continuous spatial dimensionality 
d, which has been introduced in Ref. 22, The latter is 



interesting, in particular, since in the isotropic limit non- 
trivial finite-size corrections to Tc, which differ from the 
mean-field ones, are only realized in the dimensionality 
range between three and four. In addition, analytical re- 
sultes will be given for ferromagnetic films consisting of 
a few layers. 

The main part of this paper is organized as follows. 
In Sec. H the basic equations describing the anisotropic 
spherical model and the methods of their solution are 
be briefly reviewed. In Sec. HI the analytical solutions 
for the simplified model with periodic boundary condi- 
tions (pbc) are given. In Sec. IV the available analytical 
solutions for the film with free boundary conditions are 
considered. In Sec. ^ the corresponding numerical re- 
sults are presented. In Sec. ^ the results obtained and 
some possible extensions of the model are discussed. 



II. BASIC RELATIONS 

At or above Tc in zero field, the magnetization nii 
is zero and the ASM is described by the closed system 
of equations for the correlation functions of transverse 
(a > 2) spin components, = D{maimaj) , and the 
spatially varying gap parameter, Gi. More general forms 
of these equations with m.i ^ and a lot of other de- 
tails can be found in Refs. |2^, ^ and here only 
the most important formulas will be given. In the ge- 
ometry considered, it is convenient to use the Fourier 
representation in d' = d — 1 translationally invariant di- 
mensions parallel to the surface and the site representa- 
tion in the dth dimension. For the model with nearest- 
neighbor (nn) interactions, the equation for the Fourier- 
transformed CF arm' (q) then takes the form of a system 
of the second-order finite-difference equations in the set 
of layers n = 1, 2, 00, 



26„CT„ 



{2d0/f])S„ 



where b„ is given by 



&„ = l + d[(7,G„)-i-l]+d'(l-A') 



(2.1) 



(2.2) 



and A' for the d-dimensional hypercubic lattice reads 



1 

d' ^ 

i=l 



cos(gj). 



(2.3) 
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In Eq. (pj]), 9 is the reduced temperature defined by 



T 



T, 



MFA 



(oo) 



2dJ 

IT 



(2.4) 



The quantities cro,n' and crjv +i,T t' in the nonexisting lay- 
ers, which enter equations (2.1) at the film boundaries 
n = 1 and n — N , are set to 



CTo.n' = CTAT+l.n' = 0, (fbc) 
CTo.n' — CTAT^n', 0-Ar+i^„/ = cri_„/ , (pbc) 



(2.5) 



as the free and periodic boundary conditions. The auto- 
correlation functions in each of the N layers, s„„, satisfy 
the set of constraint equations 
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70'«n(q) = 1, 



which are the consequence of the spin rigidity, |mi| — 1. 
A straightforward algorithm for the numerical solving the 
equations above is, for a given set of Gn, to compute 
all CT„„ from the system of linear equations (2.1) and 
then to put the results in Eq. ( |2.6D to obtain, after the 
integration over the Brillouin zone, the set of nonlinear 
equations for Gn- After G„ have been found, one can 
com pute the longitudinal CF crfj^, from Eqs. (2.1) and 
(2.2), where 77 is set to 1. 

For thick films far from the boundaries, it is convenient 
to consider the deviation form the bulk value 



Gl,; 



Gr, 



G< 1, 



(2.7) 



that is proportional [and for T < Tc{oo) equal] to the 
inhomogeneous part of the reduced energy pro site Uin 
(see Ref. ^). The bulk gap parameter G satisfies the 
equation 



0GP{r]G) = 1, 



where 



P{X) 



(2.8) 



(2.9) 



is the lattice Green function. The quantity Ak = 

Jk/(2dJ) for the nn interaction is given by Eq. (2.3) 
with d' ^ d and q ^ k. The solution G of Eq. (2.J) 
increases with lowering temperature 6*; at G = 1 the gap 
in the longitudinal CF closes, longitudinal susceptibility 
diverges, and the phase transition occurs. This defines 
the bulk transition temperatureO 



0,(00) = 1/P(,7), 



(2.10) 



that generalizes the well known result for the spherical 
model = 1/P(1) (Ref. |l|). The lattice Green function 
P{X) satisfies P(0) — 1 and has a singularity at X ^ 1, 
the form of which in different dimensions can be found 
in Ref. |2|. For d < 2, the Watson integral W = P(l) 
goes to infinity; thus formula ( ^.10| ) yields nonzero values 
of the Gurie temperature only for the anisotropic model. 



ry < 1. It should be noted that in the anisotropic case the 
critical indices of the model coincide with the mean-field 
ones due to the suppression of the singularity of P{r]G) 
for G ^ 1. Below 9c, the spontaneous magnetization 
appears, and G sticks to 1. 

For n ^ 1, q <ti 1, and k ^ 1, where k is the inverse 
transverse correlation length defined by 

= 2d[l/{riG) - 1] ^ 2d[l - 7/G] < 1, (2.11) 

one can reduce the second-order finite-difference equation 
Th to the differential equation for the Green's function 



dn^ 



2dGin anw = -2d9S{n - n'), (2.12) 



(2.6) where n is considered as a continuous variable, q 



\J k} + q^, and Gi„ is defined by Eq. (2/7). Note that 



the transverse bulk correlation length ^c_l = Cca = 1/k 
increases without diverging with decreasing temperature 
down to 9(. and remains constant below 0,, in accordance 
with the behavior of G. The critical-point value 



k\ = 2d[l/?7- 1] 



(2.13) 



measures the anisotropy and varies between for the 
isotropic model and oo for th e classical Ising model. Ana- 
lytical solution of Eq. (2.12) is only possible in particular 
cases, when Gi„ has a simple form and can be guessed. 
In the domain- wall geometry, Gi„ (x 1/ cosh^[(n — no)/<5], 
where no corresponds to the center of the domain wall 
and (5|-is the self-consistently determined domain- wall 
width.Lj Accordingly, the solution for Onn' is a combi- 
nation of hyperbolic and exponential functions. For the 
isotropic semi-infinite ferromagnet at the ordinary phase 



transitiouJ-p 
of Gi„ isBE 



he dimensionality range 2 < d < 4 the form 



Gin — 



d- 3 



2dn2 



(2.14) 



and Gnn' can be expressed through the modified Bessel 

functions Ifi{qn) and K^{qn). 

Note that long- wavelength equation ( ^.12 ) is valid for 
continuous lattice dimensionalities d and its solution in 
the asymptot ic re gion n 3> 1 is universal. The disadvan- 
tage of Eq. ( 2.12| ), however, is that it cannot be solved 
numerically in confined geometries, which is desirable 
in most cases (e.g., for a semi- infinite ferromagnet away 
from criticality) when there is no analytical solution. The 
problem is the stro ng divergence of Gi„ near the bound- 
aries [see Eq. ( 2.14 )] and the ensuing loss of the boundary 
condition. In fact, however, this divergence is an artifact 
of the continuous approximation. The lattice-based cal- 
culation of Ref. |2^ for the simple cubic lattice at isotropic 
criticality yields, on the contrary, rather small values of 
Gi„ near the boundary: Gn = 0.09615, G12 = 0.01254, 
Gi3 = 0.00532, Gi4 = 0.00291, etc. To get access to 
continuous dimensionalities while perserving the discrete 
structure of the lattice in the direction perpendicular to 
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the surface, a special kind of a "lattice" has been intro- 
duced in Ref. |2^. This lattice is characterized by the 
Fourier-transformed exchange interaction Jk = JoAk of 
the form 



Ak = 3 cos 

a 



^A' 



(2.15) 



[cf. Eq. (2.3)] with a continuous d' = d — 1 and in 
the whole Brillouin zone. The natural hypercubic cut- 
off < TT and the corresponding density of states are 
modified for the q components according to 



d'^\ 



i2n) 



d/_ 



dqq 



d'-l 



(2.16) 



with A = y/2{d+ 1). This choice of A preserves the 
usual sum rules for the exchange interaction. Numerical 
calculations of Ref. ^ on the continuous-dimension lat- 
tice confirmed formula ( ^.14| ) in the asymptotic region, 
n ^ 1. Near the boundary, deviations from universal- 
ity have been detected. In particular, for the continu- 
ous dimensionality d = 3.0 one obtains Gu = 0.10672, 
Gi2 = 0.01195, Gi3 = 0.00495, Gu = 0.00270, etc., 
which deviates from the results for the simple cubic lat- 
tice quoted above. 

If the gap parameter Gn is uniform, as in the bulk or in 
a film with periodic boundary conditions, Eq. (2.1) can 



be solved analytically in a straightforward way. In other 
cases, an efficient numerical algorithm is needed. One 
possibility, which is used in Ref. ^ and in the present ar- 
ticle, is to represent the layer-layer autocorrelation func- 
tion cr„„ in the continued-fraction formEjO 



2de 



1 



77 2b„ 



(2.17) 



where &„ is given by Eq. (2.2) and the functions q;„ and 
are found from the forward and backward recurrence 
relations 



1 



1 



Ctn+l = 



2K, 



26„ - a' 



(2.18) 



with boundary conditions ai — a'j^ = 0. The nondiago- 
nal CFs can then be written as 



nn' — l I I 
i=n n >n 



n 



n +1 
l—n 



ai, n' < n. 



(2.19) 



In the bulk t he qu antity b of Eq. ( p.2D is independent of 
n, and Eqs. (2.18) can be solved to give 



^ 6 - v/fo^ - 1 = 



(2.20) 



The solution the spin CFs above thus simplifies to the 
result 



O^nn' (q) = 



(2.21) 



that can also be obtained with other methods. For the 
numerical solution in the film geometry, diagonal CFs 
CT„„ of Eqs. ( 2.17 ) should be substituted into the con- 
traint equations (2.6), which defines the variation of Gn 
after application of some iterative scheme (e.g., Newton 
method). 

The re i s another routine for solving system of equa- 
tions (|2.l| ) and (2^), which avoids numerical calcula,tion 
of the integral s ove r the Brillouin zone at each step.cJ The 
solution Eq. (2.1) for diagonal CFs cr„„ can be written 
in the form 



2de 
V 



(2.22) 



wher e A ^ is the inverse of the tridiagonal matrix A of 
Eq. (|2.l|) : Ann = "^bn and An^n±i = ^1- It is convenient 
to represent matrix A in the form 



A = B 



2d'X^I, 



(2.23) 



where / is a unit matrix and _B is a tridiagonal matrix: 



Bn 



2d/ir,Gn 



B 



1. 



Then Eq. (2.22) can be transformed to 

N tt2 



2de 



2d' a: 



(2.24) 



(2.25) 



where U is the real unitary matrix (U 



Upn — Unp), that diagonalizes 13, and Ap are the eigen- 
values of B. The nice feature of this expression is that U 
and Ap are independe nt of q. Thus the integral over the 
Brillouin zone of Eq. (2.25) can be expressed through the 
lattice Green function of the layers, Pd', and constraint 
equation (2.6) takes the form 



N 



^E^P.'(2d7Ap) 



1. 



V 



(2.26) 



p=i 



Since Pd'{X) can be calculated analytically or tabulated, 
this method should be faster than the continued-fraction 
formalism described above, especially for hypercubic lat- 
tices in high dimensions. However, application of the lat- 
ter to the semi-infinite ferromagnet in Ref. ^ has shown 
that it is fast enough, the main problem being the sta- 
bility of iterations for low dimensionalities d. For this 
reason, here the continued-fraction formalism is chosen 
for numerical calculations, whereas the diagonalization 
method will be applied to find analytical solutions for 
films consisting of small number of layers. 

Before proceeding to the films with free boundary con- 
ditions, let us consider the simplified problem with peri- 
odic boundary conditions. Here, in an artificial way, all 
quantities are rendered uniform, but, nevertheless, there 
is a dimensional crossover discussed above. 
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III. PERIODIC BOUNDARY CONDITIONS 



For t he p roblem with pbc, one has G„ ~ G and 6„ = b 



in Eq. 
of Eq. 



(2^ 



) . Solving the latter with boundary conditions 



2.5) yields 



1 



(3.1) 



where a < 1 is defined by Eq. (2.20). In the long- 
wavelength region near criticality, one has 



= 1- 



(3.2) 



Thus, if NyJ + < 1, an additional factor \/ + q 



appears in the denominator of Eq. (3.1), which is respon- 
sible for the dimensional crossover in films. 

The Curie temperature of a film can be found from the 
condition that the longitudinal susceptibility diverges. 
This amounts to vanishing of the determinant of the sys- 
tem of equations for the longitudinal CF at zero wave 
vector, i.e., Eq. (2.1) with rj — 1 and q = 0. It can be 
seen that the latter is satisfied for b{ri = l,q = 0) = 1 
and hence G = 1, as in the bulk. Now 9c can be found 
from the constraint relation ( |2.6| ). Separating the bulk 
term, one can write 



d!_ 



dqq 



2d 



(3.3) 



lattice 
and 



where, for continuous dimensionalities d, the 
Green function Pdii]) is determined by Eqs. (2.9) 
( ^.16 ). For the hypercubic lattice, a usual integration 
over the Brillouin zone is performed in Eq. ( |3.3[ ). The 



large- TV asymptotes of Eq. (3.3) can be calculated ana- 
lytically. For the weakly anisotropic model in the range 
d > 3 and Nkc ^ 1, one obtains 



d' 2dT{d-2)({d-2) 



(3.4) 



For the hypercubic lattice, one should make the replace- 
ment 



d'/A"^' Sd'/{2TTy 



(3.5) 



where Sd = 27r''/2/r(d/2) is the surface of the d- 
dimensional unit sphere. Dependence of such a type, 
derived field-theoretically, can be found in Ref. For 
3 < d < 4, it is in accord with finite-size-scaling formula 
( |l . l| ) with \ ~ l/vi, = d — 2. In higher dimensions for 
the D = oo model one has l/i'b — 2, whereas the shift 
exponent A is g ivep-by he same formula. A particular 
case of Eq. (O) iMB 



Piil) + 2/{3N^) 



(3.6) 



f or d = 4. For KcN > 1, the second term in the rhs of Eq. 
(|3.3D, and thus the shift of Tc, decays exponentially with 



KcN. 



In three dimensions, the integral in Eq. ( |3.3| ) can be 
done analytically for all values oijMi^N, provided iV ^ 1. 
For d = 3 the result has the forma 



■In- 



1 



nN 1 - cxp(-KciV) 



For KcN ^ 1 it simplifies to 



P^iv) 



TtN KcN 



(3.7) 



(3.8) 



For the continuous-dimension lattice d = 3.0, the second 
terrms of the above equations should, according to Eq. 
(U), be multiplied by T[{d+\)/2\[2TT/{d+l)]^'^-^^''^ 
tt/2. The form of Eq. (3.8) is in accord with finite- 



size-scaling formula (1.1), but now the coefficient A in 



Eq. (1.1) depends on anisotropy in a crucial way. In the 
isotropic limit, Kc defined by Eq. ( 2.13| ) goes to zero, and 
the logarithmically divergent second term of this formula 
becomes dominant. In this limit 9c becomes logarith- 
mically small, in accordance with general formula ( |l.2| ), 
which for the ASM takes the form 



9-\d,r^,N)^ 



d ( dig -I 

{d-l)N '^~^\d-l 



(3.9) 



Indeed, the Curie temperature of the square lattice is 
given by 



P2(r;) = i In- 

TT 1 



^ 2 4^/2 
— = — In . 

■q TT Kc 



(3.10) 



The extra factor 3/ (27V) in Eq. (^), relative to this 
result, is due to the rescaling of the exchange interaction 
and the number of nearest neighbors in the mean-field 
transition temperature of Eq^ (^.4|). Appearance of N 



under the logarithm in Eq. (3.8) is due to cutting of the 
integral over q in Eq. ( |3.3| ) at q ~ 1/A^ rather than at 
the Brillouin-zone boundary: Only in this range of q the 

film behaves two dimensionally. 

For d < 3 the integral in Eq. (3.3) for KcN ^ 1 is 
dominated hy Kc ^ q 1/-^, and the result has the 
form 



Oc^ = Pdiv) 



d' dK' 



d-2 



h'^' KcN sin[(3-d)7r/2]' 



(3.11) 



the 



Note that because of the sma ll fac tor k'^ 
merator, both terms of Eq. ( 3.11 ) can be comparable 



with each other. For th e square-lattice stripe d — 2, the 
bulk term of Eq. ( 3.1l|) can be neglected and the result 
simplifies to 9^^ — 2/{kcN). This is in absolute accord 
with Eq. (3.9), since for the chain Pi{ri) = — i]'^. 

For the continuous-dimension model d = 2.0, the the ex- 
pression for 9~^ should be multiplied by Tr/y/6 « 1.283. 
It should be noted that the phase transition in chains 
or stripes can only occur in the limit D —t oo. For any 
finite D, the long-range order is broken by longitudinal 
fiuctuations which lead to formation of domains. 
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IV. FREE BOUNDARY CONDITIONS 

In a film with fbc, there is an additional disordering 
due to boundaries that is responsible for the larger sup- 
pression of the Curie temperature of a film in comparison 
to the model with pbc. The greatest difference between 
the two models arises within the MFA, where for the pbc 
model Tc is the same as in the bulk and for the hypercu- 
bic fbc model one hasQ 



1 



1 



1 — cos ■ 



N + 1 



(4.1) 



This result also can be obtained for the ASM in the 
limit ^ 0, where Eq. (2.1) trivializes and the mean- 
field approximation becomes exact., For the pbc model 
in this limit, finite-size corrections to Tc disappear as 
exp(— Kc-/V). 

For f] 0, the situation for the fbc model becomes 
complicated, since the gap parameter G„ is nonuniform 
near the boundaries and cannot be calculated analyti- 
cally in most cases. This problem becomes, however, 
nonessential in high dimensions, d > 4, as well as in all 
dimensions for HcN ^ 1. Here the inhomogeneity of G„ 
is localized to the narrow region near the film bound- 
aries, n '--^ 1 in the first casaH3 and n ~ 1/kc <^ N in 
the second case. In the main part of the film, G„ = G 
is uniform, and its value at the critical point of the film 
can be obtained from the condition that the longitudinal 
susceptibili ty d iverges, i.e., the determinant of system of 
equations {2A) ioi t] = 1 and q = turns to zero. It is 
instructive to calculate, ins tead of the determinant, the 
whole spin CF fro m E q. (2.1) with the free boundary 
conditions of Eq. (2.5) for the uniform G, ignoring the 
narrow boundary region where G is nonuniform. The 
result for the transverse CF reads 



,fbc. 



q) = 



de 



[l-a2"][l 



2(Ar+l-n)l 



ri\/TP 



1 _ Q,2(W+1) 



(4.2) 



the longitudinal CF being given by the sam e expression 
with Tj = 1. For the bulk CF of E n. (^.2lD , as weU as 
for the CF of the pbc model, Eq. (pj]), the divergence 
corresponds to b — 1. This leads, according to Eq. ( |2.2| ) 
with T] = 1 and q = 0, to G = 1 at the phase transition 
poin t. The pole at 5 = 1 is cance lled, however, in Eq. 
(4.2), since, according to Eq. ( 2.20| ), 6=1 entails a = 1. 
The new pole in the fbc correlation function corresponds 



N+l 



to a 



a = exp 



= —1, i.e., to 



N+l' 



-{a + a 1) 



N + l 



(4.3) 



Now the value of G at t he C urie temperature of the film 
can be found from Eq. (2.2) with t] = 1 and q = 0. The 
result has the form 



G- 



1 



1 



N+l 



(4.4) 



Now 9c can be found from the constraint relation (2.t) 
for the transverse CF of Eq. ([1.2|). The integral over 



the Brillouin zone in Eq. (2.6) simplifies in both cases, 
d > 4 and KqN ^ 1. In high dimensions, the integral 
is dominated by q ^ 1, and at such q the terms with 
a to high powers beyond the boundary regions can be 
neglected. The same situation takes place for KcN ^ 1 at 
all q; in both cases the bulk CF enters under the integral 
and the the Curie temperature of the film is given by 



= GP(77G), 



(4.5) 



with G defined by Eq. (4.4). In t he Ising limit rj 
0, the mean-field result of Eq. (4T) is recovered, since 
P{0) = 1. This result is valid for all N, since in this 
limit G„ = G is uniform throughout the film and all 
calculations above are exact. If rj substantially deviates 
from zero, this approach works only for large N, where 



Eq. (4.4) simplifies to 



(4.6) 



Since the last term in this expression is a spiall pertur- 
bation, Eq. (4.5) can be written in the formE3 



1 - 



where 



I{rj) = l + r^P'{r^)/P{v). 



(4.7) 



(4.8) 



Function Idiij) shows the same qualitative behavior as 
Pd-2{v) {see, e.g., Ref. ^). In particular, in the weakly 
anisotropic case ^ 1 [see Eq. ( [2.13 )], one has 



A-d 



/(I) - Can 



1 < d < 4 
4 < d < 6. 



(4.9) 



For d > 4, /(r/) remains finite in the isotropic limit, rj — > 
1. In this dimensionality range formula ( [4.7| ) works for 
all anisotropics, as was argues above. In the marginal 
case d = 4, 1(77) diverges logarithmically for rj ^ 1. For 
the hypercubic lattice, Eq. (4.7) takes the form 



1. 



(4.10) 



If anisotropy becomes so small that KcN <1, then the 
integral over the Brillouin zone in the constraint relation 
is cut at the lower limit q ~ 1/iV instead of q ~ Kc- This 
leads to the result 



^4(1) 



IniV 



(4.11) 



where G ~ 1 cannot be calculated analytically. For the 
continuous-dimension lattice d — 4.0, one should, accord- 
ing to Eq. (3.5), insert the factor 37r/(5vT0) ~ 1.873 in 
front of the logarithms in the above equations. 

For d < 4, I {rj) and thus the finite-size correction to Tc 
given by Eq. (4.7) diverges in the isotropic limit, 77 — > 1. 
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In the dimensionality range 3 < d < 4, however, the lat- 
ter does not imply that Tc goes to zero, since the film is 
a system of dimension d! — d — 1 and it can still order 
without anisotropy for d' > 2. In fact, the situation here 
is determined by small wave vectors, q<,l/N <C 1, in 
contrast to q ~ 1 for d > 4. In this wave-vector range d'- 
dimensionality of the film manifests itself by the appear- 
ance of an additional power of q in the denominator of 
the spin CF (T„„(q). Dimensional arguments show that 
this should be the combination Ny^ + < 1. Thus 
the finite-size correction to Tr has the form 



2d 



Pdiv) 



1 

N 



l/N 



dqq' 



rf-2 



(4.12) 



which leads to the same dependence AT^ = Tdo o) - 
Tc{N) cx l/N"^'^ as for the model with pbc, Eq. (Q 



Analytical calculation of the numerical factor in this de- 
pendence seems to be very difficult or impossible, be- 
cause it requires knowing of the self-consistently deter- 
mined variation of G„ in the film. This numerical factor 
should be larger that for the model with pbc because of 
the additional disordering at the film boundaries. One 
can see that for d < 4 this contribution to the shift of 
Tc is larger than the contribution ATc cx from the 

region q ^ I, which is described by Eq. (4.7). At large 
distances N>1/kc, the wave- vector range in which the 
film behaves d'-dimensionally disappears, as can be seen 
from Eq. (4.12). Here a crossover to the regime of Eq. 

(4.7) occurs. 

For d < 3, the intergal in Eq. (4.12) diverges at the 
lower limit in the isotropic case. In this range, the di- 
mensionality of the film is d' < 2 and it cannot order 
without anisotropy. For very small anisotropics, Tc of the 
film is very low, and at such temperatures the chains of 
spins going from one boundary to the other are strongly 
correlated and behave as single composite spins, as was 
explained in the Introduction. Here the situation sim- 
plifies, and the type of boundary conditions no longer 
plays a r ole in the determination of Tc that is given by 
Eq. (3.9). In the marginal case d = 3, the corresponding 
formula has only logarithmic accuracy,Lj and the more 
accurate result for the simple cubic lattice has the form 



^c'=P3{v) + ^ln^ 



C_ 

iV' 



(4.13) 



where C ~ 1 should be calculated numerically [cf. Eq. 
( [4.13 )] . It is interesting to note that the variation of G„ 



throughout the film can be easily found in the limit under 
consideration. It adjusts in a way that is compatible with 
the full correlation of the film's layers, i.e., with (T„„/ 
being the same for all n and n' at small q. The latter 
requires that all ctn an d in Eq. (|2.19|) are close to 1. 
Then from Eq. (2.18) one finds that 6„ = 1 inside the 
film. On the contrary, in the boundary layers one obtains 
bi = bj\i = ^. The resulting variation of Gn in the film 



can now be determined from Eq. (2.2) and has the form 



■ , , .i=l.N 

G„- <; 2d-i' 

1, n^l,N. 



(4.14) 



To conclude this section, let us consider the bilayer, 
N — 2. For the bilayer. Curie temperature can be cal- 
culated analytically with the help of the diagonalization 
formalism described at the end of Sec. ||. At Tc one has 
Gi = G2 = 2d/(2d— 1), and the Curie temperature ifself 
is given by 



E 



2d - 1 ± 



Pd-i 



2{d-l)r] 
2d - 1 ± 7? 



(4.15) 



where ^_|_ sums terms with both signs. In dimensions 
d < 3, the term with the sign minus diverges in the 
isotropic limit, rj 1, whereas the term with the sign 
plus remains finite. In this case the result is compatible 
with the general- formula (3.9). 

For films consisting of larger number of layers, the 
diagonalization formalism becomes cumbersome. For 
A = 3, in particular, the cubic characteristic equation 
for the problem simplifies, and one can write down ana- 
lytical expressions for the correlation functions in terms 
of Gi = G3 and G2. The latter and the Curie tempera- 
ture, however, are defined by a system of transcedental 
equations that cannot be solved analytically. 



V. NUMERICAL RESULTS 

The system of equations describing the anisotropic 
spherical model was solved numerically in the following 
way. For a given variation of G„, the transverse layer- 
layer autocorrelation functions (t„„ were found with the 
help of Eqs. (2.17) and (2. IS). The results were substi- 
tuted into the constraint relations ( |2.6| ) to obtain, after 
the integration over q, the system of nonlinear equations 
for the set of G„. The latter was solved by the Newton 
method. In fact, all equations were rewritten in terms 
of deviations from the bulk values (see Ref. ^), which 
helps to improve the accuracy. This calculation yield the 
variation of Gn throughout the film for any temperature 
T > Tc- To find the value of Tc, the nonlinear equation 
in temperature, I?(T) = 0, was solved, where VjT) is 

for 



the determinant of system of linear equations (2.1) 
T] = 1 and q = 0. In this sys tem of equations, the values 
of bn are given by Eq. (2^), where G„ are functions of 



temperature found on the previous step of the numerical 
routine. 

It is convenient to plot, instead of G„, its deviation 
from the bulk value Gi„ = G„ — G which at and below the 
bulk criticality (i.e., for G = 1) is equal IXLthe nonuniform 
part of the reduced energy density C/i„.c3 The variations 
of Gin in the isotropic films at bulk and film criticalities 
in four and five dimensions for the continuous-dimension 
model are shown in Fig. |l|. The results for the hypercu- 
bic lattice differ from the latter by nonuniversal factors 
of order unity and are not shown. One can see that for 
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FIG. 1. Energy-density profiles at bulk and film criticalites 
for the isotropic D = oo model in four and five dimensions. 
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FIG. 3. Reciprocal of normalized energy-density profiles at 
bulk criticalites for the isotropic D — oo model in dimensions 
d = 3, 4 and 5. 
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FIG. 2. Energy-density profiles at bulk and film criticalites 
for the isotropic D = oo model in dimension d — 3.5. 



d = 5 at the film criticality the value of Gi„ is close to 
{-k/NY /{2d) « 4 X 10"'* in the main part on the film, 
in accordance with Eq. (44). The bulk-criticality pro- 
files of Gin that are plotted for comparison, show a very 
different behavior. They follow the appropriate semi- 
infinite profiles Gin ^ Xjrfi^'^ (Ref. ^) until approxi- 
mately n ^ N/A. Crossing of the latter with the plato 
level at some n* ~ N^/i^-^) determines the boundary 
region n^n* where the semi-infinite Gin profile is real- 
ized. For d > 4, the value of n* increases slowlier that 
N] thus for thick films the plato of Gi„ dominates the 
film's Curie temperature, as was explained in Sec. 
the marginal case d = 4, the plato of Gin is mucE 
pronounced. In fact, it arises only for In ^ 1. 

For d < 4, the problem becomes more complicated, 
and no analytical solution for the variation of Gin could 
be obtained. Numerical results for d = 3.5 that are 
shown in Fig. ^ are also different for film and bulk 
criticalities. At the film criticality, the power-law de- 
pendence Gin ~ 0.165A^~^/^n^^/^ is seen, instead of a 
plato. For comparison, the dependence of Gin in the 



[V|In 

less 



semi-infinite geometry, which is given by Eq. (2.14), has 
the form Gi„ « 0.0268n~^ for d = 3.5. It is clear that for 
very thick films both criticalities are very close to each 
other, and, at the film criticality, the semi-infinite solu- 
tion should be realized in some region not too far from 
the boundary. The crossover from this solution to the 
numerically found n~'^/^ law can be located by equating 
the expressions above. The result is that for d = 3.5 at 
the film criticality, the semi-infinite solution is realized 
in the region 1 ^ n<n* with n* « 0.02647V, whereas 
the n~^^^ law holds for n* <n but not too close to the 
center of the film (see Fig. ||). What is the nature of 
the numerical parameter 0.0264? Is it small enough to 
allow an analytical derivation of the n~^/^ law? These 
questions have not been addressed in this article. 

For d < 3, the film's Curie temperature and the energy- 
density profile at the film criticality depend on anisotropy 
in an essential way. This makes the corresponding fig- 
ures less interesting; in particular, for KcN ^ 1 there 
is a plato of Gi„ at the distances from the boundaries 
n>l/Kc. By contrast, energy-density profiles at the 
bulk criticality can be studied for d > 2 for isotropic 
models. This problem for the semi-infinite and film ge- 
ometries can be addressed with field-theoretical meth- 
ods for systems with an arbitrary number of spin compo- 
nents D (see, e.g., Refs. p7| and 28). In Refs. 



29 and 3C 



the energy-density profiles in critical films with various 
boundary conditions have been calculated with the help 
of the e expansion. For the free boundary conditions con- 
sidered here, the exponentiated energy-density profile g 
is given by Eq. (5.7) of Ref. |3^. In the limit D ^ oo one 
can use v = 1/2 and v — 1 m four and three dimensions, 
respectively, and for d = 4 and d = 3 the result of Ref. 
simplifies to 



Sm TTX 



smTTX 6 



(5.1) 
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where z is the distance from the left boundary, L is the 
film thickness, x = z/L = n/N, ip{x) = r'(2;)/r(a;), 
and 7 = 0.57721566... The quantity g{x), by defini- 
tion, may differ from Gi„ only by a numerical factor 
which is dependent on N. Thus the reduced quantity 
g{x) = g{x)/ g{0.5) can be compared with the Gi„ pro- 
files numerically calculated for the isotropic D = oo 
model, which are shown in the scaled form in Fig. |^. 
One can see that for d = 3 the agreement is reason- 
ably good, keeping in mind that (73 (x) of Eq. (5T) has 
been obtained in the first order in e = 4 — d. The curve 
1/54(0;), however, goes in the middle between the curves 
representing the numerical results for d = 4 and d — 3 . 
The origin of this discrepancy is that 54(2;) of Eq. (5.1) 



[T,M-T,(N)]/U^) 



does not comprise logarithmic terms which appear in four 
dimensions. In particular_Jpr the semi-infinite problem 
in four dimensions one has£3 Gi„ = [16n^ ln(an)]~^ with 
a ^ 1 for ri ^ 1. By contrast, both 34 and g^ simplify 
to g{x) — l/x^ for x <C 1, which corresponds to the the 
semi-infinite case. For d > 4, the variation of Gi„ can 
be found analytically, at least at or above the bulk crit- 
icality, by the integration of the free proparator, as was 
done for the semi-infinite problem in Sec. IIIC of Ref. |2|. 
This routine coincides with that used in Ref. ^ in the 
film geometry up to numerical factors; thus the results 
of both approaches, in the scaled form, are identical in 
high dimensions. For d — 4, however, logarithmic term 
appears in the method of Ref. in contrast to that of 
Ref. 1^. 

Another possible normalization of the energy-density 
profiles is that with respect to the solution of the 
semi-infinite problem G^. The function g{n/N) = 
Gin/G!j"^y2 contains more information than g = 
Gi„/Gijv/2 used above. In particular, the value of g 
in the middle of the film can be used to compare the 
accuracy of different approaches. So, for d = 3 the nu- 
merical calculation yields 33(0.5) « 1.822. For d = 4, 
the value of 5(0. 5) slightly increases with TV because 
of the logarithmic corrections to scaling: One obtains 
54(0.5) « 1.772 for N = 100 and 54(0.5) w 1.795 for 
N = 200. On the other hand, the results o f R ef. |30| are 
scaled as g{x) = g{x)lA\ thus from Eq. (5.1) one ob- 
tains 53(0.5) K, 1.766 and 54(0.5) « 1.645. One can see 
that the discrepancy between the exact solution for the 
D = 00 model presented here and the field-theoretical 
results of Ref. ^ is comparable to the difference between 
the results for d = 3 and d — A. 

It is interesting to note that the numerical solution for 
d = 3 can be represented with a pretty good accuracy by 
the empirical formula 



g{n/N) 



Glr, 



(^/2)4/3 



G? 



(5.2) 



'i,N/2 sin4/3[(V2)(2x)3/2]- 

In Fig. m the corresponding curve would go exactly over 
the circles for d = 3, so it is not shown. The value 
5(0.5) = (7r/2)"'/3 ~ 1.826 is very close to the numeri- 
cally calculated value 1.822 cited above. 

Numerical results for the relative shift of Tc in isotropic 
films with free boundary conditions in high dimensions 
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FIG. 4. Curie-temperature shift in films with fbc for the 
isotropic D — oo model in four and five dimensions (continu- 
ous-dimension model). 
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FIG. 5. Curie-temperature shift in films with fbc for the 
isotropic D = oo model in dimension d = 3.5. 



are shown in Fig. |[ For d = 5 , th e data approach at 
large N the asymptotic formula (4/7) with 15(1) « 1.66 . 
For d = 4, the results for thick films follow Eq. (4.11) 
with G w 1.5. Both of these asymptotes are shown by 
dashed lines in Fig. |^. For both hypercubic lattices, the 
data are in accord with those of Ref. Rn, where a spherical 



model with separate spin constraints in each of the layers 
was used. This modified spherical model is thus better 
than the standard one using the global spin constraint. 
However, it cannot be generalized for the anisotropic case 
and applied to films in three dimensions. 

To illustrate the situation in the range 3 < d < 4 where 
the isotropic D = 00 model shows nontrivial power law 
for the Tc shift in films, the data for d = 3.5 are shown in 
Fig. p. One can see that both models with pbc and fbc 
show the same shift exponent A = d — 2. For the model 
with pbc, the amplitude A in Eq. (1.1) calculated from 
Eq. ( |3!^ ) is A « 1.788. For the fbc model, one obtains 
yl w 4 from the fit to the numerical results. This value is 
substantially larger than for the model with pbc because 
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FIG. 6. Curie-temperature shift in films with fbc for the 
ASM in three dimensions (continuous-dimension model). The 
crossover values of A'^ corresponding to KcA — 1 are shown 
by short vertical lines. 
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FIG. 7. Curie-temperature shift in films with fbc for 
the ASM in different dimensions vs anisotropy (continu- 
ous-dimension model). High- and low-anisotropy asymptotes 
are shown by dashed lines. 



of the boundary effects that affect the whole sample (see 
Fig. Also due to the boundary effects, the data ap- 
proach the power-law asymptote much slowlicr than for 
the pbc model. 

Numerical results for the Tc shift in films with fbc in 
there dimensions for different values of t he an isotropy are 
shown in Fig. |. The data follow Eq. ( |4.13| ) at n^N < 1 
and Eq. (4.7) for larger thicknesses, Kr.Nft^- In the 



latter case, the shift exponent in Eq. (1.1) is A = 2, as 
in the mean-field case, but the amplitude A can be very 
large for small anisotropy. 

Finally, the dependences of the Tc shift in films of 
A'^ = 50 layers in different dimensions are shown as func- 
tion of the anisotropy in Fig. Fo r d — 5, the data 
are well described by the formula ( |4.7| ) anisotropy range. 
For d — 4, this formula fails at small anisotripies be- 
cause of the logarithmic divergence of 14(1]) a.t rj = 1 



and it should be replaced by Eq. ( 4.11 ) in the isotropic 
limit. For d = 3.5, the Newton algorithm for finding 
Tc fails at small anisotropies because the inverse lon- 
gitudinal susceptibility of a fil m or the determinant of 
the linear system of equations (2.1) approach zero very 
fiat: Xz^{T) oc V{T) cx (T - Tc)^''"^) the isotropic 
case. Knowing this dependence, one can go over to 
the equation for the nonlinearly scaled quantity, e.g., 
jy{d-3)/2^ that behaves linearly near Tc- In such a way, 
in fact, the data in Fig. ^ have been obtained. This 
yields Tc{N)/Tc{oo) « 0.989 for d = 3.5 in the isotropic 
limit. The functional form of the dependence of Tc on 
anisotropy near the isotropic limit can be found from the 
observation that the film behaves as a d'-dimensional sys- 
tem i n the wave- vector range qN < 1. One can write [cf. 



Eq. (12.101) ] 



Tc{N) ^ Tc{N) 



Tc{oo) T'c(oo) 



1 



N 



(5.3) 



where Cd should be determined numerically. Fit to the 
data for d = 3.5 yields C3.5 ~ 2.40. With this choice, 
asymptote overlaps with the numerical results in 

Fig. |7| in a sufficiently wide range of Kc- 

In three dimensions, Tc{N) goes to zer o log arithmically 
in the isotropic limit according to Eq. ( 4.13 ). The fit to 
the numerical data for the continuous- dime nsion lattice 
[3/tt in front of the logarithm in Eq. (4.13) is replaced 
by 3/2] yields C — 1.77. Again, this asymptote well fits 
the numerical data in a wide range of anisotropies. 

For d < 3, the Curie temperature of the film goes 
to zero as Tc oc Nkc~'^ in the isotropic limit, in accor- 
dance with Eq. (3.11). The corresponding asymptotes 
for d = 2.5 and d = 2.0 are shown by dashed lines in 
Fig. ^ Unfortunately, numerical calculations could not 
been performed for such small anisotropies because of the 
instability of the numerical algorithm. 



VI. DISCUSSION 

In this paper a profound influence of the anisotropy 
on ordering in ferromagnetic films has been investigated. 
The main result is that in three dimensions, the be- 
havior of the film of arbitrary number of monolayers N 
can be dominated, for sufficiently small anisotropy, by 
two-dimensional fiuctuations which strongly suppress the 
ordering temperature. Films described by the isotropic 
Heisenberg model cannot order in three dimensions. Most 
of results above have been obtained for the exactly solv- 
able infinite-component spin-vector model, but more re- 
alistic models as the Heisenberg model should share the 
same qualitative behavior. In part icular, the small- 
anisotropy formula for Tc, Eq. (1-2), should be valid 
for all models in dimensions less than three. In this di- 
mensionality range, the shift amplitude A in Eq. (1.1) 
should diverge in the isotropic limit for all systems. 



The theoretical tool used here, the anisotropic spheri- 
cal model, allows to separate the major effects of Gold- 
stone or would be Goldstone modes in systems of reduced 
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dimensionality from usually less significant but complica- 
tive effects of critical fluctuation coupling. The ASM has 
a lot of advantages in comparison to the mean-field ap- 
proximation or to the standard spherical model, and it 
should be tried for classsical spin systems the next after 
the MFA. Some more subtle effects, as the Berezinsky- 
Kosterlitz-Thouless phase transition, are, however, not 
reproducible within the ASM. 

A peculiar feature of the ASM is that only fluctuations 
transverse with respect to the easy axis give contribution 
to the equations describing static properties of the sys- 
tem. Thus these equations contain the transverse corre- 
lation length ^c_L = 1/k as the length scale. In particular, 
crossover between different regimes for the Curie temper- 
ature of the fil m, T ^N), occurs at KcN ~ 1, where Kc is 
given by Eq. (2.12). The longitudinal correlation func- 
tion and longitudinal correlation length are the "slave" 
quantities that are decoupled from the system of equa- 
tions for the ferromagnet and can be determined after the 
solution of the latter. By contrast, in the conventional 
theory of phase transitions the diverging longitudinal cor- 
relation length is used as the length scale. It would be 
interesting to see how does it appear in the theory of or- 
dering in films, if one goes beyond the D = oo model, 
i.e., takes into account the 1/D corrections. 

A more urgent task, however, is to solve the system 
of equations describing the anisotropic spherical model 
below Tc, both in the semi-infinite and film geometries. 

To conclude, it should be noted that equations resem- 
bling those describing the ASM can be obtained for quan- 
tum systems in spirit of the random-phase approximation 
by decoupling of the high-order Green functions. Exam- 
ples of such works in the film geometry are Refs. ^ and 
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